Testing Bell's inequality with two-level atoms via population spectroscopy 
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We propose a feasible experimental scheme, employing methods of population spectroscopy with 
two-level atoms, for a test of Bell's inequality for massive particles. The correlation function measured 
in this scheme is the joint atomic Q function. An inequality imposed by local realism is violated by 
any entangled state of a pair of atoms. 
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The violation of Bell's inequality [|lj is a mathemati- 
cal expression of the remarkable quantum phenomenon 
of entanglement. It may be viewed as a manifestation 
of the irreconcilability of quantum mechanics and "lo- 
cal realism," as was first revealed by Einstein, Podolsky, 
and Rosen (EPR) in their famous gedanken experiment 

. Perhaps the most familiar version of the EPR exper- 
iment, with two spin-half particles in the singlet state, is 
due to Bohm §. 

Interest in Bell's inequalities and the EPR "paradox" 
has remained amazingly intense during many decades 
The violation of Bell's inequalities by entangled 
quantum states has been verified in a number of exper- 
iments with photons. Most of them used pairs of corre- 
lated photons emitted in atomic cascade transitions 
or produced in parametric down-conversion A seri- 
ous difficulty in experiments with photons is the limited 
detection efficiency. Nevertheless, these experiments un- 
equivocally confirmed quantum-mechanical predictions in 
a clear contradiction with inequalities imposed by local 
realism. A promising experimental scheme was recently 
proposed for testing Bell's inequality with entan- 
gled two-mode states of the quantized light field. In 
this scheme [|l^, correlation functions are given by the 
joint two-mode Q function and the Wigner function. Two 
decades ago an experiment was made [|ll[ with pairs of 
spin-correlated protons produced in low-energy proton- 
proton scattering. To the best of our knowledge, this 
was the only attempt to demonstrate experimentally the 
violation of Bell's inequality for massive particles. 

In the present Letter we propose a feasible experimen- 
tal scheme for testing Bell's inequality with entangled 
states of a pair of two-level atoms, in a situation which is 
physically very similar and mathematically equivalent to 
the familiar Bohm version of the EPR gedanken experi- 
ment. We use the fact that the interaction of two-level 
atoms with classical radiation fields is mathematically 
equivalent to the precession of spin-half particles in clas- 
sical magnetic fields. These interactions, accompanied by 
the measurement of the atomic-level population, consti- 
tute the essence of the Ramsey population spectroscopy 
method A rearrangement of this convenient exper- 
imental method was recently proposed for the measure- 
ment of the spin Q function ||l3| , in line with the ideas of 
Refs. The motivation of works was to use 



the measured Q function for the reconstruction of the ini- 
tial quantum state. Here, we employ the joint Q function 
of two atoms as the correlation function which violates 
Bell's inequality for entangled atomic states. 

If a collection of N two-level atoms (or atomic ions) 
interacts with classical light fields, one can describe this 
physical situation as the interaction of N spin-half parti- 
cles with classical magnetic fields. Denoting by the 
"spin" of the rth particle, one can define the collec- 
tive spin operators J = X^r^r- The three opera- 
tors {Jx, Jy, Jz} constitute the su(2) simple Lie algebra, 
[Jk,Ji] = i^kimJm- The Casimir operator is a con- 
stant times the unit operator, = j[j -\- 1)1, for any 
unitary irreducible representation of the SU(2) group; so 
the representations are labeled by the single index j. The 
representation Hilbert space Tij is spanned by the basis 
{|j,m)} (m = Jz\j,m) ^ m\j,m). The 
"cooperative number" j is a non-negative integer or half- 
integer and takes values ■^N, ■^N — 1, ^TV — 2, . . .. For two 
atoms, the possible values are j = 1 (symmetric triplet 
states) and j = (the antisymmetric singlet state). We 
assume that the atoms (ions) are far enough apart so their 
wave functions do not overlap and the direct dipole-dipole 
coupling or other direct interactions between the atoms 
may be neglected. 

In our discussion we will employ concepts of phase 
space and coherent states. For SU(2), the phase space is 
the unit sphere §^ = SU(2)/U(1), and each coherent state 
]l7| , p^ is characterized by a point on the sphere, i.e., by a 
unit vector of the form n = (sin0 cos (/), sin0 sin </>, cos 0). 
Specifically, the coherent states \j;n) are given by the 
action of the group elements 



5(n) 



(1) 



on the highest-weight state: |j;n) = g{n)\j,j). A useful 
phase-space distribution function is the Q function: 



Q(n) = (j;n|p|j;n), 



(2) 



where p is the density matrix of the system. 

In our scheme, a pair of two-level atoms, prepared ini- 
tially in the state 1-0), is spatially separated jl^. All sub- 
sequent operations (transformations and measurements) 
are made on each atom separately. The basic idea of the 
method can be described as follows. Each atom interacts 
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with a classical radiation field of frequency close to the 
atomic transition frequency and/or evolves freely. This 
results in a phase-space displacement (rotation) of the 
atomic state. Then one measures the population of the 
upper state 1+). Of course, such a measurement implies 
that in fact an ensemble of identically prepared systems 
is used. This scheme, which is based on the method rou- 
tinely employed in Ramsey population spectroscopy , 
gives values of the Q function for each atom. Coinci- 
dence measurements for the two atoms give the joint Q 
function. This function contains information about quan- 
tum correlations between the two atoms and can be used 
in Bell's inequality to distinguish between quantum me- 
chanics and local realism. 

In the spin description, the magnetic moment fj. = figS 
is associated with each particle. If a uniform external 
magnetic field Bq = BqZ is applied, the Hamiltonian for 
each particle is given by Hq = — /x • Bq = huJoSz, where 
hujQ — —fXoBo is the separation in energy between the 
two levels. (The two eigenstates of Hq, denoted by |-|-) 
and |— ), are the two atomic levels). The corresponding 
Heisenberg equation is dS/dt = uiqxS, where (jJq = ojoz. 
This spin precession describes the free evolution of the 
atom. In addition, one can apply the so-called clock ra- 
diation which is a classical field of the form 

Ti± — B± {y cosLut — xsinujt) , (3) 

where a; ~ and we assume cvq > 0. In the refer- 
ence frame that rotates at frequency uj, the spin S in- 
teracts with the effective field B = B||Z + B±y, where 
-B|| = Bo{ll!o — uj)/lO{). In the rotating frame, the Hamil- 
tonian is H = — /iqS • B, and the Heisenberg equation for 
S is dS/dt ^ u)' X S, where a;' = {luq — uj)z + uj±y and 
uj± = —iiQB±/h. We do not use any special notation for 
the rotating frame because the population of the upper 
state |-|-), which will be measured, is the same in both 
frames. 

In the Ramsey method there are two kinds of interac- 
tions. In the first kind, B± is nonzero and constant during 
the time interval Tg. During this period B ~ B±y, where 
it is assumed that ::§> \B^^\, i.e., \uj±\ ^ |cjo — uj\- 

(This assumption means that the free evolution can be 
neglected because it is very slow compared to the preces- 
sion induced by the clock radiation). Therefore, in the 
rotating frame, S rotates around the y axis by the angle 
9 = u±Tg. The second kind is just the free evolution. 
During the time interval the clock field B± is zero, so 
B = and S rotates around the z axis by the angle 
(f)= {ujo- uj)T^. 

Now, assume that the rth atom (r = 1,2) first evolves 
freely, resulting in a rotation around the z axis by —4>r, 
and then interacts with the clock radiation, resulting in 
a rotation around the y axis by —Or- The overall trans- 
formation performed on the whole system of two atoms 
is the phase-space displacement (rotation): 



|V'(ni,n2))= n e'''-^-e'*'-^-|V')=5l(ni)5|(n2)|^). 

r=l,2 

The last step of the experimental procedure is measure- 
ment of the upper- level populations of the atoms. This 
measurement can be made with almost unit efficiency us- 
ing the technique of "quantum jumps" pc| ]. The whole 
procedure should be repeated with many identically pre- 
pared systems. Define g^'' (?" = 1, 2) which is 1 if the rth 
atom is found in the upper state |-|-) in the Ith experi- 
ment, and otherwise. Repeating the described proce- 
dure with L identical systems (L ^ 1), one determines 
the probability to find the rth atom in the upper state, 
Qr — L^^ Y^]^=i 1r^- For the displaced state |V'(ni, n2)), 
the probabilities Qr are 

Qr{nr) = ,(+|Trp^^{|i/'(ni,n2))(V(ni,n2)|}|-f)r 

= r{nr\pr\nr}r, r=l,2. (4) 

Here, pr = Tr^^^ {|^)('0|} is the reduced density matrix 
of the rth atom, and |n,.),. — gr{nr)\+)r is the atom's 
coherent state, 

|n) = cos{9/2)e-''^/^\+) + sin{9/2)e'^/^\-). (5) 

Hence, Qr(nr) is just the Q function for the rth atom. 

In a similar way, the probability to find both atoms in 
the upper state, determined in experiments with L identi- 
cal systems, is Q12 = 1i^l2^ ■ For the displaced 
state |'(/'(ni, n2)), this probability is 

(3i2(ni,n2) = |i(+| 2(+|V'(ni,n2))|^ 

^|i(ni|2(n2|V')|'. (6) 

This is just the joint Q function for the system of two 
atoms. Therefore, "rotations" performed on the atoms 
followed by the measurement of the upper-state popula- 
tions (procedures routinely employed in population spec- 
troscopy) allow one to measure values of the Q functions 
at any phase-space point. Note that the free evolution 
(represented by a rotation around the z axis) does not 
affect the level populations, if it is not followed by a ro- 
tation around an axis different from z. 

The dichotomic observables qi and (72 and the corre- 
sponding joint probability Q12 (the correlation function) 
are just what one needs to formulate a Bell inequality. 
Define the combination of the Clauser-Horne type |2l[] 

r = Qi2(0, 0) + Qi2(n, 0) + Qi2(0, n') 

-Qi2(n,n')-Qi(0)-Q2(0). (7) 

Local realism requires that T must satisfy the inequality 
— 1 < r < 0. If two atoms are prepared in an entangled 
state, quantum mechanics predicts that a measurement 
performed on one of them "affects" the other. This con- 
tradicts local realism and is expressed mathematically by 
the violation of Bell's inequality. Consider two atoms 
prepared in the entangled state 
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This is the singlet state \ j — Q,m = 0) for = tt and the 
triplet state \j = l,ni = 0) for (p = 0. The calculation 
of the individual and joint atomic Q functions is rather 
straightforward. For 6 — 9' , we find 

r = sin^^-lsin^0cos^ ^~t'~^ -l. (9) 
2 2 2 ^ ' 

The maximum value of F is 0, so we look for a minimum. 
For (/)- (j}' = if, we get F = -1 - i(cos 9 - cos^ 9). This 
function reaches its minimum F = —9/8 = —1.125 for 
9 — 7r/3. This clearly violates the limit imposed by local 
realism. 

If one considers two light modes a and b which share 
a single photon, the description is very similar to the 
case of two spin-half particles. The singlet state then 
takes the form |u(7r)) = 2'^/^ (|l)a|0)b - |0)a|l)b), where 
the one-photon state |1) is identified with the upper spin 
state |-|-) and the vacuum state |0) with the lower spin 
state |— ). Peres |2^] proposed the following experiment: 
one observer measures the projections of |w(7r)) on the 
states |l)a and \uj+)a — ^(■\/3|l)a -I- |0)a) of the mode a 
and the other observer measures the projections of |u(7r)) 
on the states |l)b and |a;_)f, = i(v/3|l)b — |0)b) of the 
mode b. Using these measurements, one can construct 
the Clauscr-Horne combination F which violates Bell's 
inequality, F — —9/8. While in the context of the Fock 
spaces of light modes the states \uj+)a and \uj-)b look 
somewhat contrived, in the context of two-level atoms 
they are easily recognized as the atomic coherent states 
|n)i and |n')2 with 9 = 9' = tt/3 and — </>' = tt. As 
we just found above, projections on these coherent states 
can be used to obtain the maximum violation of local re- 
alism, F = —9/8. It is important to emphasize that in the 
context of two-level atoms these projections on coherent 
states are physically meaningful (and experimentally fea- 
sible) , as they can be measured by means of appropriate 
operations of population spectroscopy. 

Another possible entangled state of two atoms is 

\v{^))=2-'/m+),\+h + e^'^\-),\-h). (10) 

This state belongs to the triplet subspace {j = 1). For 
61 = 61', we find 

F = 1 (^cos9- cos^ 9 - sin^ gcos^ + ^ ^^^^ 

The minimum value of F is —1, so we look for a maxi- 
mum. For (j) + (j)' — (f = TT , we get F = i(cos 9 — cos^ 9). 
This function reaches its maximum F = 1/8 = 0.125 for 
9 = tt/3. Once again, an entangled state violates the 
limit imposed by local realism. 

It should be noted that enormous theoretical and ex- 
perimental progress has been made during the last few 
years in the generation of entangled atomic states. The 



use of trapped atomic ions was proposed in a number of 
recent works to generate the state ( |To| ) and its multi- 
particle generahzations, 2--^/^(|j, j) -f e*'^|j, — j)). Deter- 
ministic entanglement of two trapped ions was reported 
in Entangled states created in this experiment are 
similar to those of Eq. (||). Another method for the de- 
terministic generation of entangled states of two trapped 
ions was recently proposed in Experimental prepa- 
ration of entangled states of pairs of Rydberg atoms in a 
cavity was reported in ||2^ . A method was also proposed 
to entangle neutral atoms via cold controlled colli- 
sions. Maximally entangled states of three spins were 
created in a recent experiment using nuclear mag- 
netic resonance. 

It is known that factorizable quantum states (i.e., 
states with no entanglement between the two subsystems) 
satisfy Bell's inequality. Consider a system of N two- level 
atoms with "cooperative number" j which is split into two 
subsystems with cooperative numbers ji and j2 such that 
j = ji+j2 (e.g., separation of two atoms in a triplet-basis 
state with j = 1). Then the highest state and the lowest 
state factorize: \j,±j) = |ji, ±ji)i |j2, ±^2)2- More gen- 
erally, if the whole system is in the coherent state \j, n), 
then it also factorizes [ pgf : 

b,n) = g{n)\jj) = 5i(n)52(n)|ji, ji)i|i2, ^2)2 

= |ji,n)i|j2,n)2. (12) 

Here, we used the rule for the addition of angular mo- 
menta, J = Ji + J2, to obtain ^(n) = gi (11)32(11). As a 
result, the correlation function factorizes as well: 

Qi2(ni,n2)=( I ( 1 . (13) 

Therefore, a system of two-level atoms prepared in the 
coherent state never violates Bell's inequality (provided 
that j = ji + j2)- Moreover, it can be proved |2^ that 
atomic coherent states are the only pure states which fac- 
torize upon splitting and thereby never violate Bell's in- 
equality jso) . (Note that the highest state and the lowest 
state are particular cases of coherent states.) 

On the other hand, any pure entangled state (i.e., any 
pure state which does not factorize) can always, by a suit- 
able choice of local measurement apparatus, be made to 
violate Bell's inequality |^,^. In definition (|^) the mea- 
surement apparatus is partially fixed, because the mea- 
sured observables were chosen to be the upper-level pop- 
ulations of the atoms. However, we still have to choose 
displacement parameters (rotation angles), and this free- 
dom is very essential. In fact, any pure entangled state 
of two atoms violates the inequality — 1 < F < 0. This is 
intuitively clear because of the following argument. The 
measurement of the upper-level population of an atom 
gives the projection of the associated spin S on the z axis. 
However, the projection on any axis can be measured in 
this way if one forestalls the population detection with 
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an appropriate rotation. This means that in fact we are 
free in the choice of the dichotomic variables measured 
by local apparatus. 

For a more rigorous proof, let us consider the state 
|?7) = cosi5|+)i|+)2 + sini5e*'^|-)i|-)2. An arbitrary 
pure state |^) can be written in the Schmidt decomposi- 
tion as = cosi9|^+)i|C+)2 +sini?e^'^|C-)i|C-)2, where 
|^±)i and |C±)2 are some orthonormal bases for the first 
and the second atom, respectively. The state is entan- 
gled if ^ iriTT (n = 0,±1,...). Clearly, can be 
obtained from |?7) (with the same d and Lp) by rotations 
of the single-atom bases: |?/') = gi{n)g2{'ai)\rf) (the de- 
gree of entanglement, dependent on is conserved under 
these local transformations). Therefore, the joint Q func- 
tion of an arbitrary state |V') at a point (ni,n2) is equal 
to the joint Q function of |r?) at some other point. We 
are able to show that the state \r]) violates the inequality 
— 1 < r < as long as it is entangled (the amount of 
violation depends on the degree of entanglement and has 
a maximum F = 0.125 for "d = 7r/4). Consequently, any 
pure entangled state violates this inequality. 

In conclusion, we presented an efficient experimental 
method to test the violation of Bell's inequality by entan- 
gled states of two-level atoms. Our method employs the 
convenient technique of population spectroscopy which is 
characterized by very high measurement accuracy (the 
detection efficiency is almost unity). The correlation 
function measured by this method is just the joint atomic 
Q function. The corresponding Clauser-Horne combina- 
tion F violates the inequality — 1 < F < 0, imposed by 
local realism, for any entangled state of a pair of atoms. 
Maximal violations are found to be F = —1.125 and 
F = 0.125. 
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